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1 .  INTRODUCTION 


The  impact  of  potential  energy  derivative  studies  within  such  areas  as 
molecular  structure  determination,  vibrational  spectroscopy,  molecular 
dynamics,  and  statistical  rate  theories  is  only  beginning  to  be  realized. 
Statistical  and/or  dynamical  models^~^  based  upon  ideas  such  as  the  Reaction 
Path  Hamiltonian^  require  in  the  simplest  application,  second  derivatives  of 
the  potential  energy  and  in  slightly  more  refined  versions,  they  require  some 
limited  fourth  derivative  information.^  The  ultimate  impact  of  these  models 
depends  to  some  extent  on  the  availability  of  energy  derivatives  of  various 
order  for  reliable  ab-initio  wave functions. 


First  derivatives  of  the  Born-Oppenheimer  potential  energy  for  Hartree 
Fock  wavefunctions  have  been  available  for  quite  srae  time  now,  and  second 
derivatives  have  been  available  for  a  few  years.”  '  However,  in  order  to 
investigate  the  broad  regions  of  the  potential  energy  surface  implicit  in  some 
of  the  above  applications,  one  needs  a  potential  function  which  includes  at 
least  differentially,  the  effects  of  electron  correlation.  There  has  been  a 
great  deal  of  recent  activity  in  potential  energy  derivative  evaluation  for 
correlated  wavefunctions.  piere  now  exist  codes  which  calculate  analytically 
MCSCF  second  derivatives,°~^^  Cl  first  derivatives  based  upon  a  Hartree  Fock 
reference  function^^  or  on  a  multiconfigurational  reference  function^®  and 
Cl  second  derivatives  based  upon  a  Hartree  Fock  reference  function.^  Third 
derivative  equations  have  been  published  for  various  wavefunctions^^  and 
one  can  expect  to  see  the  implementation  of  some  of  these  equations  in  the  not 
too  distant  future. 


He  address  here  an  issue  pertinent  to  the  efficient  evaluation  of  any  of 
these  potential  energy  derivatives.  The  derivatives  of  the  energy  are  most 
tractably  evaluated  in  a  basis  of  cartesian  displacements  of  the  atoms  in  a 
molecule.  There  are  3N  such  displacements  where  N  is  the  number  of  atoms. 
Since  the  energy  is  invariant  to  overall  translations  and  rotations  of  the 
molecule,  the  3N  elements  of,  for  example,  the  cartesian  gradient  vector  are 
not  all  independent.  Given  3N-6  elements  of  the  cartesian  gradient  vector,  we 
show  how  one  can  generate  the  remaining  six  elements  by  consideration  of  the 
rotational  and  translational  invariance  of  the  energy.*  He  also  show  how  the 
same  ideas  can  be  applied  to  the  evaluation  of  second  and  third  cartesian 
derivatives  yielding  even  larger  relative  savings. 


II.  DEVELOPMENT 

He  denote  the  cartesian  coordinates  of  a  molecule  by  the  vector  x.  If  r^ 
is  a  coordinate  representing  an  overall  rotation  or  translation  of  the 
molecule  (an  "external''  coordinate),  then  the  change  in  cartesian  coordinates 
accompanying  a  displacement  along  r^  can  be  expressed  as 


*Fot  the  sake  of  simplicity,  the  discussion  in  the  text  Implicitly  assumes  a 
nonlinear  configuration  of  the  molecule.  The  arguments  are  fully  applicable 
for  a  linear  aiolecule  for  %ihich  one  singles  out  five  cartesian  directions 
instead  of  six. 
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Here  x  (r.>0)  are  the  cartesian  coordinates  of  the  reference  geometry. 

If  r^  represents  a  translation  of  the  molecule,  then  the  second  and  higher 
order  terms  are  aero  and  Eq.  (1)  expresses  a  linear  transformation.  If  r^  is 
a  rotation,  all  terms  contribute  to  the  sum. 

The  energy  is  invariant  to  external  displacements. 

E  (r.)  ■  constant.'  (2) 

Differentiating  Eq.  (2)  with  respect  to  r>  and  evaluating  the  result  at  r^^O 
yields  directly  a  relationship  among  the  components  of  the  cartesian  gradient 
vector. 
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where  g  is  the  gradient  vector  in  cartesian  coordinates.  There  are  six  such 
equations  corresponding  to  the  six  external  directions  r^,  i  ■  1,6.  Taken 
together,  they  constitute  a  linear  system  of  equations  of  order  six  for  which 
the  unknowns  are  six  components  of  the  cartesian  gradient  vector. 


where 
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Equation  (4)  is  the  working  equation  for  gradient  calculations.  Care  must  be 
taken  such  that  the  six  cartesian  directions  separated  out  in  Eq.  (4)  are 
linearly  independent,  otherwise  the  matrix  of  coefficients  is  singular  and  the 
solution  undefined.* 


A  useful  relationship  similar  to  Eq.  (4)  can  be  obtained  for  each  row  or 
column  of  the  cartesian  second  derivative  matrix  by  differentiating  Eq.  (3) 
with  respect  to  a  cartesian  displacement. 


*A  similar  situation  is  encountered  when  one  exploits  the  inter-relationships 
among  the  two  electron  integrals  under  a  rotation.  See  L.R.  Kahn, 
"Relationships  Among  Derivatives  of  the  Integrals  in  the  Calculation  of  the 
Gradient  of  the  Electronic  Energy  with  Respect  to  the  Nuclear  Coordinates,"  J. 
Chem.  Phys..  Vol.  75,  p.  3962,  1981. 
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In  this  equation,  F^  is  the  k^^  row  or  column  of  the  cartesian  aecond 
derivative  matrix  and 

tJJ  -  (L  )  .  (7) 
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There  are  six  of  these  equations  just  as  in  the  gradient  case.  Casting  this 
equation  in  the  same  working  form  as  Eq.  (4),  the  vector  F  takes  the  place 
of  g  and  the  right  hand  side  now  has  added  to  it  the  second  (gradient 
dependent)  term. 
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It  is  interesting  to  note  that  L  is  an  eigenvector  of  F,  with  eigenvalue 
zero  only  if  the  second  term  in  Eq.  (6)  is  zero.  This  is  always  true  for 
translations  since  in  this  case  T  -  0.  For  rotations,  this  is  the  case  if 
the  gradient  is  zero,  i.e.,  at  a  stationary  point. 

The  extension  to  third  derivatives  is  straightforward.  Eq.  (3)  is 

differentiated  with  respect  to  two  cartesian  components  X|^,  x^  to  obtain 

relationships  among  the  elements  of  the  vector  G, 

^kl 


where 


G.  ,  -  ^  (F. -) 
~kl  dx  kl 


One  obtains 


^  -.1  -It  -  -It ,  1 

L  G.  ,  +  T  F,  +  T  F,  +  T  g 

~r.  ~kl  ~r.  ~k  ~r.  ~1  ~r.  » 

1  111 


•li  L.  G...  -  -  b’^^  i  -  1,( 
J“1  J,r.  jkl  1 


-•  .  •  A 


(12) 


-  ^5_  L.  G.  ♦  (T^  F  ♦  F  +  g.)  (13) 

1  J*7  j,r^  jkl  j«i  j,r.  jk  j,r.  jl  J 


T*'-'  ■  ,  (L  ) 

-r.  JX|^3Xj  ~r. 


The  linear  equations,  Eq.  (12)  can  then  be  solved  for  each  k,  1  pair. 


III.  IMPLEMENTATION 

As  previously  mentioned,  the  six  cartesian  directions  which  are  separated 
out  must  not  be  linearly  dependent.  To  insure  this,  we  choose  them  as  all 
three  components  of  one  atom,  two  components  of  a  second  atom  and  one 
component  of  a  third  atom  which  is  not  collinear  with  the  first  two  atoms. 

The  three  vectors  L  corresponding  to  translations  can  be  chosen  as 

. 

normalised  translational  vectors  in  the  three  cartesian  directions.  The  other 
three  vectors  can  be  chosen  as  rotational  vectors  about  any  three  linearly 
independent  axes.  For  convenience,  we  choose  these  to  be  the  cartesian 
axes.  Thus  to  first  order,  a  rotation  about  the  s  axis  for  example  is  a 
cartesian  vector  for  which  the  y  component  on  each  atom  is  proportional  to  the 
X  coordinate  of  that  atom  and  the  x  component  is  proportional  to  the  negative 
of  the  y  coordinate.  Three  such  vectors  are  formed  and  subsequently  Schmidt 
orthogonal ized  to  the  translations.  The  result  is  three  first  order  rotations 
in  the  cartesian  planes.  Infinitesimal  displacements  along  the  direction  of 
these  vectors  correspond  to  infinitesimal  rotations  of  the  molecule. 

The  T  is  the  derivative  of  the  external  vector  L  with  respect  to  a 
~r.  ~r. 

cartesian  displacement  and  can  be  readily  generated  by  a^finite  difference 
procedure . 
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The  T  *  can  be  generated  similarly. 

~r . 

1 

IV.  DISCUSSION 

Using  the  procedures  outlined  above,  one  can  ignore  six  of  the  cartesian 
directions  in  derivative  evaluations.  These  components  are  calculated  from 
the  others  by  s  simple  procedure  at  the  end  of  the  calculation.  It  is 
worthwhile  to  examine,  for  a  few  specific  cases,  just  what  computation  is  by 
passed  by  using  this  technique. 
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The  most  straight  forward  derivatives  to  evaluate  are  those  based  upon 
wavefunctions  which  are  determined  variationally .  Recognizing  the  stationary 
condition  with  respect  to  the  variational  parameters,  one  can  avoid  the 
derivatives  of  variational  parameters  with  respect  to  nuclear  coordinates. 

The  difficult  step  in  SCF  and  MCSCF  first  derivative  calculations  arises 
through  the  dependence  of  the  energy  on  the  positions  of  the  basis  functions; 
the  evaluation  of  the  derivatives  of  the  two  electron  integrals  with  respect 
to  nuclear  coordinates.  Using  the  technique  presented  here,  one  can  select 
six  directions  for  which  derivatives  of  the  integrals  need  not  be 
calculated.  Since  the  translational  invariance  of  the  integrals  themselves  is 
exploited  in  most  derivative  integral  computer  programs, the  actual  savings 
one  can  expect  comes  from  the  neglect  of  three  cartesian  directions.*  One  can 
choose  these  components  to  be  on  the  atoms  having  the  largest  number  of  basis 
functions  and  therefore  affect  the  largest  number  of  integrals. 

Gradients  for  Cl  wavefunctions  based  upon  SCF  or  MCSCF  orbitals  and 
second  or  third  derivative  for  SCF  or  MCSCF  wavefunctions  require  the  first 
order  response  of  the  SCF  or  MCSCF  variational  parameters  to  nuclear 
displacements.  This  leads,  for  each  cartesian  displacement,  to  the  evaluation 
of  a  linear  system  of  equations  the  order  of  which  is  the  number  of 
variational  parameters.  Solving  these  coupled  perturbed  Hartree-Fock  or  MCSCF 
equations  can  become  a  time  consuming  step  in  these  calculations.  Six  of  these 
linear  systems  can  be  avoided. 

Cl  second  derivatives  and  MCSCF  or  SCF  fourth  derivatives  require  the 
second  order  changes  in  the  SCF  or  MCSCF  variational  parameters.  This  leads 
to  a  linear  system  of  equations  for  each  pair  of  cartesian  coordinates.  One 
need  solve  only  1/2  (3N-6)  (3N-5)  instead  of  1/2  (3N)  (3N+1)  sets  of  linear 
equations.  Cl  second  derivatives  have  the  additional  time  consuming  step  of 
solving  the  coupled  perturbed  -  Cl  equations.  This  requires  an  integral 
transformation  and  solving  a  linear  set  of  equations,  the  order  of  which  is 
the  number  of  Cl  coefficients.  Thus,  the  solution  of  the  coupled  perturbed 
-Cl  equations  for  each  cartesian  coordinate  is  approximately  equivalent  to  the 
evaluation  of  the  original  Cl  wavefunction.  For  relatively  small  systems, 
avoiding  six  of  these  coordinates  yields  a  significant  reduction  in  computer 
time. 

The  emphasis  in  this  discussion  has  been  on  the  application  of  the 
described  technique  to  analytical  derivative  evaluations.  This  is  because 
there  is  a  great  deal  of  current  activity  in  this  area.  One  can  also  apply 
this  scheme  to  derivatives  calculated  by  finite  difference  procedures. 


*Rotational  relationships  among  the  two  electron  integrals  can  also  be  used  in 
integral  first  derivative  evaluation.  However,  considering  the  complexity  of 
the  Implementation  in  conjunction  with  the  number  of  operations  required  to 
solve  the  equations  versus  solving  the  integrals  themselves,  this  is  not 
generally  done. 


In  summary,  we  have  presented  equations  expressing  the  interrelationships 
among  the  cartesian  derivatives  of  the  Born-Oppenheimer  potential  energy. 

These  relationships  are  derived  from  the  invariance  of  the  energy  to  overall 
rotations  and  translations  of  a  molecule.  We  show  how  the  energy  invariance 
can  be  used  to  reduce  the  computational  effort  for  various  derivative 
calculations.  Finally  it  is  emphasized  that  these  equations  are  both  simple 
to  implement  and  themselves  require  negligible  computational  labor. 
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